%load("wavelet.mat")
CA1 =CA(:,1:16);
CA2 =CA(:,17:24);
CA3 =CA(:,26:29);
CA4 = CA(:,29:30);
CA5 = CA(:,31);
CD1 = CDD(:,1:16);
CD2 = CDD(:,17:24);
CD3 = CDD(:,26:29);
CD4 = CDD(:,29:30);
CD5 = CDD(:,31);
s=zeros(1,1);
for i = 1:1795
% [CA1(i,:), CD1(i,:)] = dwt(s, 'db1'); % 采用db1小波基分解
% a1 = upcoef('a', ca1, 'db1', 1, Len); % 从系数得到细节信号
% d1 = upcoef('d', cd1, 'db1', 1, Len); % 从系数得到细节信号
% s1 = a1+d1; % 重构信号
% [CA2(i,:), CD2(i,:)] = dwt(CA1(i,:),'db1');
% [CA3(i,:), CD3(i,:)] = dwt(CA2(i,:),'db1');
% [CA4(i,:), CD4(i,:)] = dwt(CA3(i,:),'db1');
% [CA5(i,:), CD5(i,:)] = dwt(CA4(i,:),'db1');

CA4(i,:) = idwt(CA5(i),CD5(i),'db1');
CA3(i,:) = idwt(CA4(i,:),CD4(i,:),'db1');
CA2(i,:) = idwt(CA3(i,:),CD3(i,:),'db1');
CA1(i,:) = idwt(CA2(i,:),CD2(i,:),'db1');
 temp=   idwt(CA1(i,:),CD1(i,:),'db1');
 s(i)  = temp(32);
%
% figure;
% subplot(2, 2, 1); plot(ca1); title('ca1 低频');
% subplot(2, 2, 2); plot(cd1); title('cd1 高频');
% figure;
% subplot(2, 2, 1); plot(ca2); title('ca2 低频');
% subplot(2, 2, 2); plot(cd2); title('cd2 高频');
% figure;
% subplot(2, 2, 1); plot(ca3); title('ca3 低频');
% subplot(2, 2, 2); plot(cd3); title('cd3 高频');
% figure;
% subplot(2, 2, 1); plot(ca4); title('ca4 低频');
% subplot(2, 2, 2); plot(cd4); title('cd4 高频');
% figure;
% subplot(2, 2, 1); plot(ca5); title('ca5 低频');
% subplot(2, 2, 2); plot(cd5); title('cd5 高频');
% figure;
% subplot(2, 2, 1); plot(ca6); title('ca6 低频');
% subplot(2, 2, 2); plot(cd6); title('cd6 高频');

end


